EPJ manuscript No. 

(will be inserted by the editor) 



o ! Finite Size Scaling of the Spin stiffness of the Antiferromagnetic 
^ S = \ XXZ chain 

a ■ 



Nicolas Lanorencie, Sylvain Capponi and Erik S. S0rensen 



[ Laboratoire de Physique Quantique & UMR CNRS 5626, Universite Paul Sabatier, F-31062 Toulouse, France 



c3 



c 

o 
o 

CN 
> 
CN 

^1" 

o 



i 

c 

o 
o 



X 
S3 



February 1, 2008 

Abstract. We study the finite size scaling of the spin stiffness for the one-dimensional s = ~ quantum 
antiferromagnet as a function of the anisotropy parameter A. Previous Bethe ansatz results allow a deter- 
mination of the stiffness in the thermodynamic limit. The Bethe ansatz equations for finite systems are 
solvable even in the presence of twisted boundary conditions, a fact we exploit to determine the stiffness 
exactly for finite systems allowing for a complete determination of the finite size corrections. Relating the 
stiffness to thermodynamic quantities we calculate the temperature dependence of the susceptibility and its 
finite size corrections at T — 0. A Luttinger liquid approach is used to study the finite size corrections using 
renormalization group techniques and the results are compared to the numerically exact results obtained 
using the Bethe ansatz equations. Both irrelevant and marginally irrelevant cases are considered. 

PACS. 7 5.10.-b, 75.10.Jm, 75.40.Mg 



1 Introduction 

The general XXZ Hamiltonian on a ring of size L is given 

by 



H = H + Hi(A) = J - J2(S+Sr +1 + he) + J A £ S*S* +1 . 

(1) 



i=i 



1=1 



Here H is the free part and Hi(A) is the interacting part. 
It is well known that this model is solvable when peri- 
odic boundary conditions are applied jjj. However, the 
same model can also be solved under more general bound- 
ary conditions, in particular under the so called twisted 
boundary conditions defined by: 



S 



L + l 



C2 

°1 ; 



°L + 1 



where tp is the twist angle. The application of such bound- 
ary conditions is equivalent to considering a system threaded 
by a magnetic flux of strength ip(hc/e) [H. This fact was 
exploited by Shastry and Sutherland ||j7j] to calculate 
transport properties of the system. Notably, the total ground 
state energy as a function of ip, E (p), can be calculated 
and hence also the spin current and spin stiffness. In the 
thermodynamic limit they showed that the spin stiffness 
is given by: 



P_ 

J 



7r sin(/i) 
4 mO- p)' 



cos(/x). 



(3) 



Hence, for the antiferromagnetic Heisenberg model, A = 
1, p/J = 1/4. In the context of mesoscopic physics it is 



quite interesting to study transport properties for finite 
systems where coherence effects are important [|J and no- 
tably the finite size dependence of the current, suscepti- 
bility and stiffness are important for a complete under- 
standing of the experimental results. The stiffness, p, of 
the system is closely related to the dc conductivity of the 
system and a non-zero p implies quasi long-ranged correla- 
tions with power-law decay. In the present paper we show 
that it is possible to calculate numerically exactly the spin 
stiffness, p(L), for a finite system. p(L) can then be used 
to calculate the susceptibility for finite systems and finite 
temperatures. Using a Luttinger liquid approach it is pos- 
sible to understand quite completely the structure of the 
finite-size corrections and we compare these perturbative 
results to the numerical ones. 

We consider exclusively the antiferromagnetic (AF) 
case with J = 1 and we shall mainly be concerned with 
the regime where the anisotropy parameter, A, lies be- 
tween A = (XY) and A = 1 (XXX). It is well known 
that in the regime A E [—1,1] this model display gapless 
excitations with power law correlations and off-diagonal 
long-range order characterized by a non-zero stiffness, p, 
in close analogy to the superfluid order parameter in the 
two-dimensional classical XY model. For A > 1 the model 
Eq. (Q) enters a phase with Ising like AF order and a non- 
zero gap. Following the above remarks, this transition, 
occurring at A — 1, can be viewed as a metal-insulator 
transition Q. We concentrate on the region A S [0, 1] with 
emphasis on the finite-size corrections at A = 1. We first 
briefly review the finite-size scaling of the stiffness in sec- 
tion ||, then we discuss the simple free case, A = 0, (XY). 
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Section |] presents the renormalization group (RG) ap- 
proach leading to the predictions for the finite-size scaling 
behavior discussed in section [5] considering both irrelevant 
(A £_[0, 1[ ) and marginally irrelevant (A — 1) cases. Sec- 
tion p-R present the numerically exact method for calcu- 
lating the stiffness (and implicitly the susceptibility) and 
the numerical results up to L = 10000 are compared to 
the RG predictions from the previous sections. 



2 Scaling of Stiffness 

Suppose a twist of size tp is applied at the boundary of an 
otherwise uniform system. It is natural to expect that this 
will give rise to a uniform phase gradient 



V0 = <p/(aL), 



(4) 



throughout the system with L sites and lattice spacing 
a. We can now define the stiffness with respect to the 
resulting change in the ground-state energy density in the 
following way : 



5e 



1 



;P(V0) 2 , 



(5) 



where Seo is the change in the ground-state energy den- 
sity when the twist p is applied. It then follows that the 
stiffness is given by the following expression : 



P(L) 



(aL) 2 d 2 e Q (p) 
h dp 2 



(6) 



Hence, p has dimension of inverse (length) xt. In the 
vicinity of a quantum critical point or line we can invoke 
hyperscaling and two-scale factor universality (see refer- 
ence B for a discussion and references) to argue that 



P e~Hr = c, 



(7) 



where C is a universal constant. Applying standard finite 
size scaling theory |l(J we then expect the stiffness to obey 
the following finite size scaling ansatz M: 



p(L) = L- 



l ~p(L llv 6) 



(8) 



where S is the distance to the quantum critical point. 
In a phase with long-range order, the stiffness should di- 
verge and in the absence of long-range order we expect p 
to vanish exponentially with the system size. Hence, Eq. 
@ describes the finite-size corrections close to a critical 
point. In the present case of the one-dimensional Heisen- 
berg chain we expect to have z = 1 and consequently 
— d — z + 2 = 0. Since the twist is applied in the XY 
component of the spins we expect the resulting stiffness 
to be non-zero and universal in the critical region between 
A e [—1,1], Eq. (||) then tells us that the leading finite- 
size corrections in this region are absent and in the ther- 
modynamic limit we expect p to jump discontinuously at 
A = —1,1, as noted in Ref. ||. It is important to note 
that the above finite-size scaling analysis do not take into 
account corrections to scaling which, as we shall see later 
on, are especially important close to A = 1. 



3 Free case, A — 

At A = 0, H — Ho and the Hamiltonian ([|) becomes 
equivalent to a system of free fermions that can be di- 
agonalized explicitly in k-space through the use of the 
Jordan- Wigner transformation S? = 1/2 — rij, and 5j — 

Ej—i 
i=i ni . The ipj satisfies fermionic commutation re- 



lations, { -0| , -0j } 



Sij, and rij 



iptipj is the number of 



fermions (spin down) at the j-site. We can now write Hq 
in k-space 



H = -jy"^cos(fc)y>|y>fc. 



(9) 



The antiferromagnetic ground state is in the S* ot = sec- 
tor, which corresponds to half-filling (n to t — L/2). The 
effect of the twist tp is a shift in the momentum space: 
k — ► k + tp, this can also be seen explicitly from the 
Bethe ansatz equation (see section^). Therefore, we ob- 
tain a rather simple expression for the ground state energy 
per site versus tp : 



eo(L,p) 
J 



1 cos(g) 



(10) 



Using Eq. 
are then : 



), the finite size corrections for the stiffness 



h 



1 



L sui(jt/L) 



The corrections of order o ( -pr ) can be understood in terms 
of conformal field theory [jllj. We also have the universal 
corrections to the free energy which are given by &o{L) = 
eo° — fr§ + ®{t?)- Here, u is the velocity of excitations 
(it = 1 in the free case) and c is the conformal anomaly 
number; c = 1 for the s = 1/2 XXZ chain 



4 Renormalization group analysis 

4.1 Bosonization of the XXZ chain 

In this section, we introduce the bosonization of interact- 
ing fermionic systems like the XXZ chain (using Jordan- 
Wigner transformation) or the Hubbard model, which are 
Luttinger liquids In a few words, it consists in lin- 

earizing the dispersion relation near the Fermi points and 
studying the low energy excitations close to those points. 
Moreover, for commensurate fillings of order n, umklapp 
scattering (n left electrons become right or vice-versa) is 
possible and gives rise to an additional Sine-Gordon term 
in the Hamiltonian. In the continuum limit and in the 
Sl ot = sector (corresponding to half-filling) , the XXZ 
chain is governed by a Hamiltonian expressed in term of 
conjugate boson fields <P(x) and n(x) as follows : 



H = 



^{(^)l«H( X )f 

2 7T 



(f)[^l 2 } 
K 7T 



JA 

2{iraf 



dxcos(A$(x)). 



(12) 
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The last term represents umklapp scattering ans u, K are 
the two Luttinger liquid parameters which describe the 
low-energy physics. 

If the ring is threaded by a flux ip, the effect is just a 
shift of the current II(x) like II(x) — > LT(x) + @. 
Since we are dealing here with fermions carrying a charge, 
we can identify the factors entering the Hamiltonian as the 
charge stiffness, given by — , and the compressibility — . 
Note that we are working now in units where % = 1, a = 1. 

Using Jordan- Wigner transformation, we obtain, in 
the spin chain case the spin stiffness, 



and the susceptibility, 



uK 



K 

^ UTT 



(13) 



(14) 



The susceptibility can then be obtained from u and p using 
the hydrodynamic relation fcfl: 



P 

X = — ■ 



(15) 



In the thermodynamic limit, we will note u* and K* for 
the Luttinger liquid parameters. By using the thermody- 
namic relations from above and comparing to exact results 
obtained from Bethe ansatz, we can express these param- 
eters versus A = cos(/x) : 



J 

K*(A) = 



ir sin(/i) 
2 Ji ' 

2(tt-a*)' 



(16) 



For 1/2 < K < 1 (A £ [0, 1[ ), the perturbation g 
remains irrelevant and vanishes rapidly with the size 
L, and K goes to K* . 

At the isotropic point K = 1/2 (A = 1), the perturba- 
tion is marginally irrelevant implying a logarithmically 
decreasing g and K — > K* . 

For K < 1/2 (A > 1), the umklapp term, g, increases 
with the system size because the perturbation is rele- 
vant and K vanishes. A gap opens up in the spectrum 
and the stiffness falls to zero. 





1 s 




^^^^ K- 











figure 1. RG flow of the Sine-Gordon model. 



4.2 RG equations 

Now, we have to study the Sine-Gordon term in Eq. ( p"2| ) 
(g J dx cos(4:<P(x))) as a perturbation and examine the ef- 
fect of a renormalization under a change of length scale. 
We then obtain : 



dl 
dK 

~df 



= (2-4K)g + 0(g 2 ) 
= -Ag\ 



(17) 



where g — JA/(2(ira) 2 ) is the coupling, I = ln(L) and A is 
a constant. Those equations are identical to the Kosterlitz- 
Thouless renormalization group analysis (jl5| used in the 
classical 2d XY model and which gives a description of the 
superfluid transition. There is a powerful analogy between 
the superfluid density which vanishes at the transition and 
the spin stiffness in the quantum XXZ chain. Here, the 
renormalization is done versus the chain length L at T = 
0, but due to the invariance of the model it should be 
equivalent to considering an infinite chain at T = u/L 
|l6| . The renormalization flow diagram for this model is 
shown in Fig. [l] where distinct behaviors are observed : 



5 Finite size scaling 

The integration of the RG equations (|l7]) from Lq to L 
gives us the finite-size scaling of K(L) and g(L). The cor- 
rections to g have already been evaluated by Cardy [|l7] 
and by Lukyanov fll8| |. They are very important to cal- 
culate the corrections to the energy |12, lE,n9|. At the 



isotropic point, the coupling g is marginally irrelevant and 
the integration of Eq. ( |l7| ) up to 0(g 2 ) gives rise to loga- 
rithmic corrections to the ground-state energy |19] : 



e (L) 



^-[c + ^bg{Lf 



(18) 



where c = 1 and b — 4/V3 is determined by a three-point 
function p2]|. Here, the corrections of K will imply correc- 
tions to the stiffness and susceptibility in two cases: when 
the perturbation is irrelevant and marginally irrelevant. 



5.1 Irrelevant perturbations: power-law corrections 

In the anisotropic antiferromagnetic regime (1/2 < K* < 
1), we can integrate Eq. ( [l7|) from Lq to L, assuming the 
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integration with L 3> Lq 
corrections to K : 



We then obtain the finite size 6 Bethe Ansatz Equations 



K(L) = K* 



(K(L )-K*)(1-2K*) 



1 - K* - K(L ) 

-16(iC*-l/2) 



-8(K*-l/2) 



(19) 



We see that the exponent of 1/L lies between + and 4. 
In the free case (K = 1), we find an exponent 4 which 
seems to contradict the finite-size corrections found in 
Eq. (11). However, it must be noted that terms in 1/L 2 



the so-called 'analytical' corrections, are expected to ap- 
pear in all quantities since they are related to the confor- 
mal symmetry of the fixed-point Hamiltonian Jll],[l4|]. We 
can also emphasize that these analytical corrections dom- 
inate when K > 3/4 (i.e A < 1/2) and are sub-dominant 
otherwise. This will be shown precisely with the numerical 
Bethe ansatz calculation. 

From Eq. (|l^), we expect that to the lowest order, the 
finite size corrections for the stiffness and the susceptibil- 
ity are : 



If 1/2 < A < 1 : 

p(£)- — ~ X(L) 
If < A < 1/2 : 

,(£)- — - X(L)- 



K* 



K* 

UTT 



&{K*-l/2) 



(20) 



5.2 Marginal irrelevant perturbations: logarithmic 
corrections 

At the isotropic point, we can integrate Eq. ([l7|) at the 
lowest order, taking initial conditions at Lq. We find that 
K is decreasing logarithmically : 



K(L) = \ 



K(Lq) 



l + 4(K(L )-l)ln(±) 



(21) 



Such behavior is in agreement with the logarithmic cor- 
rections obtained by Loss et al. pOfl . At first order in 
l/ln(L/Lo), the finite size scaling correction for the stiff- 
ness and the susceptibility are 



p(L )/J- 1/4 



l + 8(p(L )/J-l/4)ln(L/L )' 



(22) 



The solution of the one-dimensional antiferromagnetic Heisen- 
berg model due to Bethe |lj can be extended to quite 
general boundary conditions. Hamer, Quispel and Batch- 
elor Q exploited this fact to calculate the ground state 
energy in the thermodynamic limit for the Hamiltonian 
with twisted boundary conditions, Eq. (^|), as a function 
of the applied twist p. Using this result the stiffness of the 
system can be calculated [|| using the relation Eq. ([}]). If 
twisted boundary conditions are imposed the expression 
for the ground state energy per site remains unchanged 
and are the same as for the uniform case: 



e (L,p) A 



L/2 

-^cos(fc i ((p)). 



(24) 



/=i 



The change in the boundary conditions enters in the equa- 
tions determining the L/2 quasi- momenta, fc/, only in the 
following manner: 



(25) 



where 0i tTl is given by: 
6>/,„ = 2 arctan[ 



Z\sin(^ 



k, -k„ 



cos(^Mp) -Z\cos( 



kl—k-n 



) 



(26) 



For the ground state, the set of integers is given by = 
I— (i/2 + l)/2, [1, L/2]. Using the above equations it 
is possible to numerically determine the quasi-momenta, 
ki, even for relatively large systems and for non-zero tp. 

Using Eq. (||) and Eq. (^4| ), we can now formally write 
an equation for p(L): 



P{L) 



L/2 



d 2 k 



J ^ L dip 2 dp 



isin(fc0 + (^) 2 cos(fc,)]| v=0 . (27) 



Hence it is possible to rather easily calculate p{L) if the 
derivatives d 2 ki/dip 2 and dki/dp> can be calculated. For 
dki / dp> we obtain the following expression: 



dp 



1 

L 



L/2 



dO(k h k n ) dki dO{k u k n ) dk n 



dh 



dp 



dk n dp 



(28) 



JX(L) 



J X (Lo) - 1/tt 2 



l + 27r 2 (J X (Lo)-lA 2 )hi(L/L )' 



(23) 



As noted previously, these results also give the tempera- 
ture dependence of these quantities by taking T = u/L 
and we obtain results in agreement with previous work by 
Eggert et al. fl6l| . These predictions will be compared to 
numerically exact results in section 0. 



Since the derivatives d0{ki,k n ) / dki do not depend on 
dki j dp, but only on the previously determined quasi-momenta 
ki, this is a simple matrix equation from which dki /dp can 
be determined using standard linear algebra routines. An 
equivalent expression exists for d 2 ki/dp 2 which also re- 
duces to a linear algebra problem once dki /dp is known. 
Hence, p(L) can be determined numerically exactly once 
the ki's have been obtained. 
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figure 2. The spin stiffness p(L) of the antiferromagnetic chain 
as a function of the anisotropy A for different ring sizes, L. The 
solid line is the result Eq. (H). In the inset is shown the uniform 
susceptibility versus A obtained from Eq. ((Fa) 



figure 3. The numerically determined exponent of the lead- 
ing correction to scaling in 1/L of p(L) as a function of the 
anisotropy A (circles). The solid line indicates the RG result 
for this exponent 8(K* - 1/2) = *Ztc^(a) ■ 



7 Numerical results 

We begin by a discussion of our results for the spin stiff- 
ness, p{L), for finite rings of size L. In Fig. our results for 
the spin stiffness are shown as a function of the anisotropy 
A for different ring sizes, L. As noted in section ||, we ex- 
pect the leading finite size correction to be absent in the 
gapless regime A £ [0, 1[. This is rather clearly the case. 
Only when A > 1 do these corrections become important, 
the system opens up a gap and the correlation length, £, 
becomes finite. In this regime, an explicit expression for 
£ has been obtained by Baxter and due to this finite 
correlation length we expect p{L) to vanish in an expo- 
nential manner with the system size, p(L) ~ e- L /Z. For 
A < 1, the finite size corrections are expected to have the 
form (^0|). Finally, at A — 1, we expect logarithmic correc- 
tions of the form (|2|) . These logarithmic corrections are 
non-negligible and remain sizable for systems of macro- 
scopic size. Hence, the discontinuous jump in p at A = 1 
is difficult to observe. 



7.1 Power-law corrections 

We now turn to a discussion of our results for the finite 
size corrections in the critical region A < 1 . Following our 
results from the previous sections, we expect p(L) to be 
independent of L to leading order and the corrections to 
scaling to acquire a power-law dependence with an expo- 
nent depending on K* . 

Writing p{L) = p + 0(L~ a ), with p given by Eq. (|), 
we can determine the exponent a from the numerically 
exact results obtained from the Bethe ansatz. We have 



determined this exponent in the region A < 1, and our 
results are shown in Fig. ||. The solid line indicates the 
renormalization group result, Eq. (|l9|). As previously ex- 
plained we expect on general grounds always to have cor- 
rections to scaling of the form 1/L 2 , as explicitly shown 
in Eq. ([TT|). Corrections of such a form dominate when 
< A < 1/2. Accordingly, the corrections to scaling com- 
ing from the marginally irrelevant coupling are dominant 
when 1/2 < A < 1. 



7.2 Logarithmic corrections 

Finally we discuss our results for the isotropic point. At 
A = 1, the perturbation to the fixed point Hamiltonian is 
marginally irrelevant, giving rise to finite size corrections 
of logarithmic form for the spin stiffness, Eq. (^), and 
the susceptibility, Eq. (|23|). In order to compare our nu- 
merical results to the RG results we use our results for the 
largest possible system size as the reference point, Lq, and 
perform the comparison for L < Lq Taking Lq = 10000, 
we show in Fig. | our results for stiffness as a function 
of system size, L. The solid line indicates the RG result 
for the stiffness, Eq. (p2|), and the numerical results from 
the Bethe ansatz solution are shown as circles. Numerical 
calculations have been performed for system sizes up to 
10000 sites. A good agreement between the RG results and 
the numerical results is evident even down to rather small 
system sizes. From Eq. (|l5| ) we see that the susceptibility 
follows a similar form. 

Using the relationship L < — > u/T, we can obtain in- 
formation about the temperature dependence of the sus- 
ceptibility in the thermodynamic limit from our results 
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figure 4. Spin stiffness, p(L) as a function of 1/L at A — 1. 
The RG result Eq.(§|) (solid line), with L = fOOOO is com- 
pared to exact numerical results from the Bethe ansatz (o). 
The thermodynamic limit value, p(l/L = 0)/J = 1/4, has 
been obtained from Eq. (^). 
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J % 0.105 
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- RG 

o Bethe Ansatz 




0.002 
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figure 5. Susceptibility x(T) as a function of the temperature 
T. The circles indicate the exact Bethe ansatz results, the solid 
line is the RG result Eq. @ with To = tt/(2 * Lo) and the 
dashed line the result from Ref. |l8| . The inset shows the tem- 
perature dependence of the susceptibility for different values of 
the anisotropy A. For A^\, the susceptibility follows a power 
law form with an exponent of 2 at A — and 0.4. At A — 0.7 
the exponent is close to 8(K* — 1/2) ~ 1.356 as expected. 



obtained at T — for finite systems. The temperature 
dependence of the Drude weight has recently been cal- 
culated using the thermodynamic Bethe ansatz |22| and 
is currently a topic of discussion [^3| . Previously, the low 
temperature behavior of the susceptibility has also been 
studied in detail by Eggert et al. jl6| and Lukyanov fl8f . 
In analogy with the work of Eggert et al. m& , we obtain 



from Eq. ( p3| ) the following expression for the low temper- 
ature dependence of the susceptibility in the thermody- 
namic limit : 



JX(T) 



J X (T ) - 1/tt 2 



^ l + 27r 2 (Jx(To)-l/7r 2 )ln(T /T) 



(29) 



As previously done, we use our largest system size to de- 
fine an equivalent temperature, T , and use the RG ex- 
pression Eq. (|2^) for higher temperatures (smaller system 
sizes). Taking T = u/L = ttJ/(2 * L Q ) ~ 1.57 x 10~ 4 
(L = 10000), we show in Fig. | the RG result, Eq. (§|) 
(solid line) along with the exact numerical results from the 
Bethe ansatz (circles) and the result of Ref. ||l8| . As al- 
ready shown by Eggert et al. p6f for the isotropic case, the 
agreement between the RG result and the Bethe ansatz re- 
sults is excellent at low temperatures (large system sizes) . 
However, the result of Lukyanov (ll| seems to work bet- 
ter over the complete temperature range. The logarithmic 
dependence of the susceptibility at A = 1 is clearly visible 
and the derivative with respect to temperature of x(T) 
diverges as T — ► and x(T) — > 1/tt 2 . 

As pointed out in the previous sections, we expect a 
power-law dependence of x(T) with the temperature for 
A < 1 with an exponent of 2 for A < 1/2 and a non-trivial 
exponent in the regime 1/2 < A < 1. In the inset of Fig. [|, 
x(T) is shown for several different values of A < 1. For 
A > 1 the system opens up a gap and the susceptibility 
decreases exponentially at low temperatures. 



7.3 Coupling term g 

From our results for the susceptibility (or equivalently 
the stiffness) it is possible to obtain an estimate of the 
flow of the Sine-Gordon coupling term g(L) with the sys- 
tem size L. Previous studies have obtained the same flow 
from the ground state energy [fl2|| . At the isotropic point 
{A = 1), we have seen that the perturbation term in 
Eq.(fl2"|), g J dxcos(4^(x)), is marginally irrelevant. Hence, 
we expect a logarithmic behavior of g(L). From the results 
of Eggert et al. JTg) for the susceptibility in the k = 1 
WZW non-linear a model 0] , the coupling g can be ex- 
pressed in term of x(L) obtained numerically from the 
solution of the Bethe ansatz equations as follows: 



BAE 



7T\/3 



(ML) 



(30) 



As mentioned, it is also possible to estimate g{L) from 
the finite size corrections to the ground state energy p2| , 
Eq. (Q). In this case one finds: 



12L 2 /7r 2 (e^-e Q (L))-l 



327T7V3 



1/3 



(31) 



where eg(L) again is determined from the numerical solu- 
tion of the Bethe ansatz equations. Finally, we can com- 
pare these two estimates to the RG result for g(L) which 
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is given by the same expression than Eq. (PQ) except that 
x{L) is given by the RG expression Eq. Q23] ) : 



J X (L ) - 1/tt 2 



2 l + 27r 2 (J X (L )-l/7r 2 )ln(L/L )' 



(32) 



In this equation, we use x(-^o) to determine the bare cou- 
pling g$. We can also compare our results to the general 
expression, including higher order corrections, derived by 
Lukyanov J18[ for g(L) which is, in our notation: 
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(33) 



where 7 is the Euler constant. In Fig. 
for gcg E {L) (solid diamonds) and g BA 



line indicates the RG result, Eq. (132 
indicates the result of Lukyanov [[18 



we show results 
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figure 6. The flow of the coupling constant g(L) with sys- 
tem size L for the isotropic Heisenberg model. The solid line 
indicates the RG result, Eq. (^). The symbols indicate esti- 
mates of g(L) from the scaling of the ground state energy and 
susceptibility and the dashed line represents the result Eq. B3. 



size, Lq — 10000, has again been used to determine the 
bare parameters in Fig. ^. As evident from the results 
shown in Fig. ^ there is an excellent agreement between 
the numerical estimates for g(L) obtained from the ground 
state energy en(i) and the susceptibility, xiLd an d the RG 
result, Eq. ((32|). The result of Lukyanov agrees very 
well with the numerical results over the entire range of L. 



A = 1 where logarithmic corrections are dominant, our 
results are in good agreement with previously known re- 
sults. From the finite size dependence of the spin stiffness 
and susceptibility we are able to obtain precise results 
for the low temperature behavior of the susceptibility in 
the thermodynamic limit. The described effects should be 
observable in mesoscopic systems where a detailed under- 
standing of the finite size behavior is crucial. 
We thank P. Azaria for useful discussions. 
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8 Conclusion 

The finite size corrections to the spin stiffness and subse- 
quently the susceptibility have been studied in detail for 
the 5 = 1/2 antiferromagnetic spin chain as a function 
of the anisotropy parameter A. At the isotropic point, 



